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Resonances in positronium hydride 

Joseph Di Rienzi, College of Notre Dame of Maryland, Baltimore, MD 21210 
Richard J. Drachman, NASA Goddard Space Flight Center, Greenbelt, MD 20771 

We re-examine the problem of calculating the positions and widths of the lowest-lying 
resonances in the Ps + H scattering system which consists of two electrons, one positron and one 
proton. The first of these resonances, for L=0, was found by the methods of complex rotation and 
stabilization, and later described as a Feshbach resonance lying close to a bound state in the closed- 
channel e + + H "system. Recently, results for the L=1 and 2 scattering states were published, and it 
was found, surprisingly, that there is a large shift in the positions of these resonances. In this work 
we repeat the analysis for L=1 and find an unexpected explanation for the shift. 

I. INTRODUCTION AND HISTORY 

The first evidence for the existence of at least one resonance in the low-energy scattering of 
positronium (Ps) on hydrogen atoms came from a stabilization calculation [1] followed by complex 
rotation computations [2], At the time it was not clear what the mechanism of the resonance might 
be, but the authors suggested that it was connected with the degenerate 2s-2p threshold in the Ps 
atom that lies nearby. Because of the neutrality of the systems this was not likely. 

A better suggestion was that there should be resonances generated by the Rydberg bound 
states in the closed, re-arranged e + + H -channel. There should be an infinite number of such states, 
and they should lie fairly close to the levels of hydrogen, although perturbed by the extended charge 
distribution of the hydrogen negative ion. We can also expect that the L degeneracy of the hydrogen 
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levels should be broken by the short-range repulsive potential of the H~ ion so as to raise the low-L 
states more than the higher ones, since the low-lying states are more penetrating. Furthermore, these 
resonances should be found only in the electronic singlet state, since the hydrogen ion is in a singlet 
state. It should be possible to carry out a fairly simple calculation based on this picture in order to 
obtain values for the resonant energies (presumably shifted somewhat from the Rydberg positions 
due to coupling with the open channel) and the widths. 

This program was carried out [3] using a simple form of scattering trial function, but one that 
included a reasonable representation of both the open and closed channels. In that calculation it was 
found that there was only a very small shift of the first two s-wave resonances away from the 
Rydberg energies, and the width of the first of these was in fairly good agreement with earlier 
calculations. This encouraging result made it seem that the mechanism of resonant scattering in the 
Ps+H system was well understood and that no surprises would be found in the higher angular 
momentum states. 

Recent results [4] of a very extensive complex-rotation calculation forL=l and 2 have turned 
out to be quite different. There is a large shift upward of the lowest-lying p-state resonance from its 
expected location near the 2p Rydberg energy, and a corresponding, but smaller, upward shift of the 
lowest-lying d-state resonance. The position of the p-state resonance is actually higher than the first 
s-wave resonance, reversing our qualitative expectation based on the picture discussed above. The 
authors of Ref. [4] suggested that this result could be explained by assuming that there is a very large 
coupling shift in these higher angular momentum states, although its cause remained unknown. 

This Paper is the result of an investigation into this puzzling situation, again using the 
methods of Ref. [3] extended to higher angular momentum. We are not trying to improve on the very 
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accurate values obtained in [4], Rather, we are trying to understand those somewhat counterintuitive 
results. First, we will formulate the scattering method used and will re-calculate the s-wave results, 
verifying the earlier work. Then the p-wave problem will be examined, and an interesting qualitative 
conclusion will be reached. The d-state problem will be relegated to a future report. 


II. FORMULATION OF THE METHOD 

To describe both the open channel [ F(R) ] in the conventional way and the closed channel 
[G(x)] in terms of bound states of the positron in the field of the hydrogen ion we write the Ps-H 
scattering function in the following form: 

Y(r p r 2 d0 =F(R,)<J)(p 1 )iKr 2 ) + F(R 2 ) <f>(p 2 ) t(^) +G{x)$(r v r 2 ). (1) 

In this expression r x and r 2 are the position vectors of the two electrons relative to the proton taken 
as fixed at the origin, and xis the position vector of the positron. The relative coordinates describing 
the position of the center of mass of the Ps atom and its internal coordinate are 
R = (r. +x)/2 and p ; =x-r, respectively.The internal wave functions of the H atom and the Ps atom 

are, of course, known exactly: 

_p 

2 

4r(r) = — and <f>(p) = — — • ® 

The nonrelativistic Hamiltonian for the interacting system (in atomic units with energies m 


Rydbergs) is 
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H=-V\-V ’l~vl+- 

x r n r, r 2 p, p 2 


( 3 ) 


It is convenient to break up the Hamiltonian in different ways corresponding to each of the channels 
in the trial function. When considering the open channel the following breakup is useful: 


where 


K, = - 2 vl r 1, Hh - -vj-l. V « = 2_ + 2 




7 r. 


r !2 X r !.2 P 2 .l 
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For the closed channel we break the Hamiltonian as follows: 


H=H , 


where H h - 




o 7 ^ 

-v;-v 2 -— 

r i 


V + W 


r, r. 


V 


12 
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* P ! P 2 
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We could treat both i 7 and G in the trial function as completely free to be determined variationally 
[5], but in the spirit of the earlier calculation we will represent G as a finite sum of bound-state 
functions: 

G(x) = C n G n^’ where / f f dr 1 d? 2 dx V 2> H > r 2 ) = € n +E H ' ' ( 6 ) 
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Since $ is the wavefinction of the negative hydrogen ion, it is an eigenfunction of (defined in 
Eq. (5)) with eigenvalue £’ H -=-1.0555, and Eq.(6) becomes 
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We can solve for these energies and wave functions both variationally in the form 


GJx) 


~-x L e~ ax 


P l (x)E d > 


W k 

A . 


£=0 


( 8 ) 


and by solving the variational differential equation numerically. 

-VX (?) + U(x) G n (x) = e n G n (x) . ' ( 9 ) 

(Later we will discuss the details of this part of the program, including approximate forms for O and 

hence U(x) and the low-lying eigenvalues e n .) 

Without loss of generality and for simplicity we will next derive the optical-potential 
scattering equation with only one term in the expansion for G. We begin with the usual functional 

I=ffJdR 1 dp l dr 2 '¥(H-E)'¥. ( 10 ) 

Using the wave function from Eq. (1) with Gix) = C n G n (x), where n refers to a particular eigenstate 
supposed to be the main component of the resonance of interest, we can write the compact form 

W=P+C Q, and I=A+BC+DC „ 2 , 
where A=(P[ff-£[P), + {P\H-P\Q), 


(ID 
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We first set the variation on the coefficient C n to zero, obtaining the result C n = -BUD . But from 
Eqs.(6) and (1 1) we see that D =E H - +e n -E ; this has the form of a resonant denominator. After having 
accounted for the coefficient C n we can write 1 in the following form: 

/-(/ y-tft. W-jWVP-W . (12) 


Varying in the usual way on the function F we obtain the following scattering equation: 


i r , . , , , V (R) fdR'v (R)F(R) 

-I[v 2 +* 2 ]F(tf) +fdR'Wji)F(R')* 2J2 -Z_ i “ 


E E H- e n 


= 0 , 


( 13 ) 


where E H = -lRy, £- ps = -^Ry, k 2 =2(E-E ps -£’ H ) = (2£'+3)Ry 


(Here we have restored the sum over bound states.")' The kernel ■X r 'comes“from“exchange“ofthed:wrf 
electrons; it is a very complicated expression and was evaluated long ago by Fraser [6], The 
potentials V n appearing in the separable integral terms are defined as follows: 


V n (R x ) = f f [H-E] G n (x)<S>(r v r 2 ) 

= ] 7^ 1 ^(P 1 )^(r 2 )[F fl -t/ + e n+ £ ff --£']G^(r 1 ,r 2 ), 


( 14 ) 


where we have used the separation of the Hamiltonian from Eq. (5) and the differential equation of 
Eq. (9) in the second form of V „ appearing in Eq. (14). Eq. (13) has exactly the form of the well- 
known Feshbach optical-potential equation [7] with one important exception: We have not made use 
of projection operators in our derivation, so the open- and closed-channel parts of the wave function 



-7- 


of Eq. (1) are not orthogonal. For that reason, the coupling potentials V n contain explicit energy 
dependence, not present in the strict Feshbach formalism, as can be seen from Eq. (14). 

In the next Section, we will discuss simplified forms of the kernel K and the hydrogen ion 
wave function <E>, and we will construct and display the coupling potentials. 


III. EXPLICIT FORM OF THE SCATTERING EQUATION 
A. The static-exchange approximation 

Consistent with our program of continuing the type of calculation introduced in [3] we will 
not try to use the exact form of the kernel K given by Fraser [6], but instead will replace it by a 
separable kernel of the following simple form (for L=0 and 1). 

K L (R,R)=g L P L (R)P L (R )e^ iR+R \ (15) 


The static-exchange approximation corresponds to the first two terms of Eq. (13), and it has been 
numerically solved [8] with the exact exchange kernel K for L=0,l,2. For each partial wave we can 
solve the corresponding scattering equation with our separable kernel analytically and will adjust the 
constants a and g to give good agreement with the results of Ref. [8]. The idea is that the adjusted 
kernel can then be used unchanged in the full optical equation. For L=0 and 1 we make the usual 
substitution F, (R) = ^ PAR) and find the following integro-differential equation to be solved: 


R 


-ul\R)-k\(R)+y^u L (R)+X L Re' aLR fdRRe' aLR u L (R), where 


8ng L 


( 16 ) 


R 


Solutions for L=0 and 1 can be easily obtained, and they have the following forms. For L-0. 



- 8 - 


For L=l: 


u Q (R ) = sin(£/?) + tan(r) 0 ) 


cos(M)-<r aoJ? 

- , \1 
1 J<x 2 +k 2 ) R 


2a. )_ 


and /tcot(rt 0 ) = 


{a Q +k 2 ) 4 k 6 +5k A a 2 Q +l5k 2 OL 0 -5a.Q 


4 Vo 


16a n 


(17) 


a,J? 


u^R) =j ! (kR) - tanfrjj) n^kR) + e 1 


1 a a 2 +k 2 „ , 
— +— + R\ 

kR k 2k 


where j (z) = S * n ^ -cos(z) and n.(z) = - C0S ^ -sinfe), 

z z 


(18) 


and ^co^rij) = - 


(a 2 +k 2 ) 4 (k 2 -a{){a*+l0a]k 2 +k 4 ) 

+ 


4X 


16a 


When we fit these analytic expressions to the numerical values given in [8] we obtain the following 
values of the constants X L and a L :X 0 =-2. 18372, a 0 =l. 10432, Aj = -2.31168, and 0^=1.08183. It 
is interesting that these constants are almost the same for both values of L. 

B. Constructing the e + -H' bound states 

Next we must choose a reasonable form for the function 0 appearing in the definition of the 

potentials V n [Eq. (14)] that represents the hydrogen negative ion. It would, of course, be best to use 

some sophisticated expression, perhaps in the Hylleraas form with many linear parameters. 

Following the simple ideas of Ref. [3], however, we continue to use the single-term, closed-shell 

wave function 3>(r. ,r 0 ) = — e Tl+r2 , which leads to simplified mathematics. This function is not 

TX 

a very accurate one for the H ~ system, since it actually does not even give rise to a bound ion, but 
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it does mimic the distribution of charge density fairly well. In any case, it is the first function that 
one should try; later other simple forms may be worth exploring. (We have already made use of the 
correct binding energy of the ion, E H '=- 1.0555, regardless of the actual variational energy 
corresponding to this form of <£.) 

The potential U in which the positron moves then has the following form [Eq. (7)]: 


U(x) = ^-ffdr l dr 2 e 

7l~ J J 


-2Z(r r r 2 )l 





( 19 ) 


(We let Z=ll/16, the best variational value.) Notice that this potential has the correct asymptotic 
form, becoming -2/x for large x values (giving the infinite set of hydrogenic states discussed above) 
and becoming +2/x for small x values (the repulsive short-range potential that should raise the low-L 
states more than the higher-L ones.) In Ref. [3] a table of the first few energies obtained by solving 
Eq. (9) are given for L=0,1, and 2. Four of those energies are of interest here:e ls =-0.3653, e 2s =- 
0.1449, e 2p =-0.2121, and e 3p = -0.0990. These correspond to Ps scattering energies (k 2 /2=0.4445+e D ) 
of 0.0792, 0.2996, 0.2324, and 0.3455, respectively. 

It has always been understood that the Is bound state should not be associated with a 
resonance but is just a contribution to the non-resonant scattering and also lowers the energy of the 
PsH bound state. The 2s scattering energy should be compared with the value 0.2944 obtained by 
Yan and Ho [9]; the agreement is quite close. But the 2p scattering energy given above is not at all 
close to the first scattering resonance energy 0.3149 given in Ref. [4], Actually, our 3p value is 


closer. 
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C. The resonant separable potentials. 

We are interested in solving Eq. (13) for p-wave scattering with the separable form of the 
kernel K given in Eq. (15) and with the potentials V 2p and V 3p calculated using the simple product 
form for the hydrogen ion wave function. By using this function in Eq. (14) we can carry out the 
integration over dr 1 analytically: 


/■ 


df n 


~(Z+l)r 2 


ypK [ 


f(x,p)-—\ 

8/jt 

P 2 J 

(i + z) 3 ; 




where f(x,p) U(x)+e +1.5 ~ . 

x Pi 2 


( 20 ) 


The remaining integral over dp { can be done by numerical integration over the magnitude of p t and 
/x = R I 'P 1 , after we re-write the positron bound-state function as follows: 


G Np (x) =g(x)x-z 



( 21 ) 


Here z refers to the fixed external coordinate system. Since both x and r 3 depend on fi, it is 
convenient to make the change of variables 


P{Z =Pi~ E Y^iR^YriR.vO-p^Rpz. 


( 22 ) 


m =- 1 


(The expression beyond the arrow can be used in our integral because the terms in the summation 
with m^O will vanish since the rest of the integrand involves only m=0.) Finally, we have two terms 
in the potential, one of which depends on the scattering energy: 


V Np (R) = A Np (k) vil ] (R)+V$(R), where A Np (k) = 0.4445 +e Np 


V, 


{1.2}. 


Np 


'N 


24 


71 


Z + l 


3 00 1 __P 

i?-f J dp j'dfie 

0 -1 


2 r + 2H\ e 


-Zr gW I 1 ] 

2 ) x |^(x)p 2 -2pj 


(24) 


and q(x) = e 


= - -(Z + D* 


f 2 


+Z+1 


+ --4e 
x 


-2Zr 


X • 


In carrying out the two numerical integrals, one must substitute the explicit forms of x and r in terms 
of R, p, and fx : 


X 

r 



2 


+— ±Rp[X. 

4 


(25) 


Then, for a set of values of R (up to the reasonable upper limit of 15) we carry out the indicated 
double integration, after which we interpolate to obtain smooth forms of the potentials for use in 
solving the separable integro-differential equation (13). 


IV. SOLUTION OF THE SCATTERING EQUATION 

The complete equation we must solve is an extension of Eq. (16), including both the term 
involving the exchange kernel K and as many “resonant terms as desired. 


dh4R) 

dR 2 


-k 2 u(R) 


2u(R) 

R 2 


-KRe + 


\6nR 


E 

N=2 


^Np 


(26) 
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In this equation the integrals I N are defined in terms of the function u(R ) as follows: 


I lp =jdRRe ~ a ' R u(R\ and l Np = f dRR V Np (R)u(R), for N>1 , 

o 


( 27 ) 


and A Np =—-e Np -0. 4445. The non-iterative method of solution of such equations works by first 
making the following decomposition of the solution: 


u(R) - % (S) + I lp u lp (R) * £ I «P U Rp m . 

-* N ^ Np 


( 28 ) 


Then the complete solution u(R) is obtained from the partial functions u Np (R) that satisfy the 
following equations: 


Qu lp (R)+X 1 Re~ aiR = 0 
Qu Np (R)+RV Np (R)= 0, [N= 2,3,...], 

where Q = -- — -k 2 +—. 

dR 2 R 2 


( 29 ) 


Two of these equations have simple analytic solutions: 


%(R) =j l (kR) and u lp (R) 


2 Jck, 


ip y 


2\2 


(ct x +k 2 ) 


n 1 ( kR ) + e 


( 2 , N 

-a,J? 1 «i *l +k 

v kR k 2k ) j 


-R 


( 30 ) 
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(These are obviously related to the solution of the static-exchange approximation.) The other 
inhomogeneous differential equations must be solved numerically. When u(R) is inserted in the 
definitions of the integrals I np a set of inhomogeneous linear algebraic equation for the quantities l np 
is obtained and is easily solved. This finally gives the scattering function explicitly, from whose 
asymptotic form the phase shift is read off. 

V. RESULTS AND DISCUSSION 

We began our investigation of the p-wave scattering by including only the lowest “resonant” 
term obtained from the 2p state. The results are shown in Fig.l, where it is clear that there is no sign 
of any resonant behavior at all. We searched quite thoroughly near the expected position of the 
resonance (k=0.6818) and also near the position found for the lowest resonance in Ref. [4] 
(k=0.7936) without finding anything but smooth behavior of the phase shift. On this figure we have 
also plotted the static-exchange results, and it is clear that adding the 2p term has raised the non 
resonant phase shifts quite considerably. As a last resort we reduced the coupling potential by 
various factors as small as 0. 1 and were able to produce a narrow, almost unshifted resonance of the 
expected kind. Our tentative conclusion at this stage was that the coupling between the closed and 
open channels had broadened the expected resonance so radically that it could no longer be 

identified; for all practical purposes the resonance had disappeared. 

We then replaced the 2p term by the 3p term and recalculated the phase shifts. These are 
shown in Fig. 2, where clear resonant behavior is seen. A least-squares fit was made to the standard 


resonant formula 
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T1 (E) 


=A+BE + arctan 

. 


r 

2 (E r -E) ’ 


( 31 ) 


and this gave an excellent fit to the computed phase shifts. The resonance we obtain 
(E r = 0.3648 Ry and r=0.0339 Ry), however, is wider and higher in energy than that calculated 
in Ref. [4]: E R = 0.3149 Ry and r=0.0032Ry. (In this case the coupling with the open channel has 
raised the resonance position by 0.0193 Ry, and its width is about 10 times as wide as in Ref. [4].) 
Finally, we included both 2p and 3p terms and obtained the results also shown in Fig. 2; the 3p 
resonance remains but is now shifted downward in energy and narrowed. The best fit was obtained 
with the following values of the resonance parameters: £ R =0.3287 Ry, and r=0.0060 Ry, which 
compare more favorably with the values reported in Ref. [4] than those determined from the 3p state 
alone. (The non-resonant parts of the 2p-3p curve are essentially unchanged from the curve with 2p 
alone.) 

Our final conclusion from this is the following: The lowest-lying p-wave Ps+H scattering 
resonance is not generated by the lowest-lying bound hydrogenic p-state in the e + +H re-arranged 
channel, as had been assumed in Ref. [3].In fact, it is the first excited p-state [3p] that produces the 
lowest resonance, while the 2p term is effective only in increasing the non-resonant scattering phase 
shift. So the apparent large shift between the closed-channel bound-state energy and the position of 
the lowest true resonance is spurious; the character of the resonance was just not correctly identified. 
(There is no reason to believe, however, that there should not be an infinite number of higher 
resonances corresponding to the infinite set of higher Rydberg states. ) 
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We have referred to this result as unexpected, as indeed it is. But perhaps (in hindsight) the 
result should not have been quite so surprising. After all, exactly the same phenomenon occurs in 
the s-wave case: the Is Rydberg state does not produce a resonance either, and the 2s term generates 
the lowest-lying s-wave resonance. In that case, the existence of the true bound PsH state may have 
confused the issue, leading us not to expect a Is resonance (In the Appendix we briefly re-examine 
the s-wave calculation.) We see now, however, that this phenomenon may be more widespread. In 
particular, we will soon be extending the present work to the L=3 case, where there is a 
corresponding, but smaller, discrepancy [4]; perhaps the same explanation will be applicable there 

as well. 
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APPENDIX 

Before beginning the p-wave calculations described in this Paper, we wished to review the 
s-wave problem that had been treated in [3]. The intent was to try to verify the results found there- 
that the coupling between the closed-channel 2s wavefunction and the open-channel scattering 
function shifted and broadened the resonance by a small amount. This would justify the reliance on 
the energies of these bound states as good approximations to the resonance energies, at least for s- 
waves, and perhaps also for higher angular momenta. It was not that we doubted the old results, but 
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we planned to use the coordinate-space integro-differential equation method described here, rather 
than the momentum-space integral-equation technique employed previously, and it seemed 
worthwhile to check the two methods against one another. 

The new results are essentially in agreement with those of Ref. [3]. In addition, by a minor 
modification of the scattering calculation we obtained the binding energy of the PsH state. With the 
inclusion of the Is and 2s terms the binding energy is increased from 0.01686 Ry in the static- 
exchange approximation to 0.03007 Ry. (This should be compared with the accurate PsH binding 
energy of 0.078393 Ry [9].) Again we find no resonance near the position of the Is closed-channel 
bound state. A least-squares fit of the present s-wave phase shift gives the following values of the 
resonant parameters: 1^=0.2955 Ry and r=0.0019, compared with the values in [3] of 1^=0.2954 
Ry and T=0.0022 Ry and the precision results in [9] of E R =0.2944 Ry and F=0.0070 Ry. The 
agreement between the various values of the resonant energy is quite satisfactory, while the resonant 
width is not as good. The main point, however, that the lowest bound state does not produce a 
resonance, is still supported by our new calculation. (As in the 2p case, we tried reducing the 
amplitude of the Is potential and recovered an almost unshifted resonance.) 
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FIGURE CAPTIONS 

Figure 1. P-wave phase shifts as functions of the Ps scattering energy. The solid line 0 ives the 
static-exchange results [7] fitted according to Eq. (18). The points (+) are the results when only 
the 2p closed-channel function is included; there is no resonant behavior seen. 

Figure 2. Resonant behavior of the P-wave phase shifts as functions of the Ps scattering energy. 
The solid points (■) are the results when only the 3p closed-channel function is included. The 
other points (+) give the results when both 2p and 3p functions are included; the resonance is 
lowered and narrowed. The curves are least-squares fits to the analytic form given in Eq. (31). 







